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CN ! Abstract 

^/y . Through multiple-scales and symmetry arguments we derive a model set of am- 

Ph ! plitude equations describing the interaction of two steady-state pattern-forming 

P3 ' instabilities, in the case that the wavelengths of the instabilities are nearly in the 

■ ratio 1 : 2. In the case of exact 1 : 2 resonance the amplitude equations are ODEs; 
I ^ i ! here they are PDEs. 

I . We discuss the stability of spatially-periodic solutions to long-wavelength distur- 

^ I bances. By including these modulational effects we are able to explore the rele- 

■ vance of the exact 1 : 2 results to spatially-extended physical systems for param- 
(<~^ . eter values near to this codimension-two bifurcation point. These new instabilities 
T-H I can be described in terms of reduced 'normal form' PDEs near various secondary 

' codimension-two points. The robust heteroclinic cycle in the ODEs is destabilised 

. by long-wavelength perturbations and a stable periodic orbit is generated that lies 

i close to the cycle. An analytic expression giving the approximate period of this orbit 

•'— I ' is derived. 

a : 
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1 Introduction 



Pattern-forming instabilities occur in many physical problems, for example 
Rayleigh-Benard convection, Faraday wave experiments and directional solid- 
ification [8]. In some of these situations well-established governing equations 
are available which are sufficiently simple to analyse. In others the situation 
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is not so clear cut, and reductions to model equations are of great value. The 
derivation of model equations to describe instabilities in these physical prob- 
lems often provides a clear and unified viewpoint, bringing out similarities in 
the underlying mathematical structure. 

In one spatially-extended dimension, the study of the stability of spatially- 
periodic patterns that arise from a steady-state bifurcation with continuous 
translational symmetry often reduces to the investigation of an evolution equa- 
tion, for example the Ginzburg Landau equation. Such equations, derived by 
asymptotic methods rather than rigorous analysis (though see recent work by 
Melbourne [17]) can both comprehend uniform spatially-periodic patterns and 
describe their stability to long-wavelength disturbances. 

In the present paper we consider the case when two distinct instability mech- 
anisms are present, for example in two-layer thermal convection [23]. Here the 
two instabilities will typically have different preferred horizontal lengthscales 
given by the critical wavenumbers corresponding to local minima in the curves 
of marginal stability. Each instability separately can be described by a 'uni- 
versal' model equation such as the Ginzburg-Landau equation. However, the 
two instability mechanisms interact nonlinearly if the instabilities occur for 
similar values of the control parameters; this is a codimension-2 bifurcation. 
In such a situation more complicated model equations are needed to describe 
the dynamics in the region of parameter space near the codimension-2 point. 
The form of these model equations will depend on the symmetry of the prob- 
lem, the steady or oscillatory nature of the instabilities, and on the ratio of 
the critical wavenumbers. 

When the critical wavenumber ratio is rational (say p/q where p and q are 
coprime integers) we may restrict attention to solutions which are periodic in 
the horizontal direction and derive, in a mathematically rigorous fashion us- 
ing a centre manifold reduction, ODEs describing the dynamics on the centre 
manifold at the point of instability. The cases p = 1, g = 1,2,3 have been 
termed 'strong' spatial resonances since the nonlinear interaction terms gen- 
erated appear at third order, or below, in the resulting amplitude equations. 

When the ratio is irrational this centre-manifold approach cannot be applied. 
Any restriction to spatially-periodic solutions will not be able to capture im- 
portant features of the dynamics. In consequence the periodic solutions found 
may be unstable to long-wavelength instabilities, and their subsequent evo- 
lution must be described by PDEs rather than ODEs. In this case, there is 
currently no rigorous mathematical formulation which leads to PDEs. At best 
the equations used represent asymptotic approximations to the true situation. 
Nonetheless we shall adopt the PDE approach here, following the work of 
(among others) CouUet & Repaux [5], since it is natural and has proved very 
productive in similar, though simpler, problems. 
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In this paper wc discuss the interaction of two steady-state instabihties with 
wavenumbers close to, but not exactly in, the ratio 1 : 2. using multiple- 
scales expansions in time and space. The case of exact 1 : 2 resonance was 
treated first by Dangelmayr [9] and later by Jones Sz Proctor [13] and by 
Proctor & Jones [23] (hereafter referred to as Part 1) in the context of two- 
layer thermal convection. Important results were also obtained by Armbruster 
et al. [1] and Julien [12], and more recently by Porter & Knobloch [20]. In 
particular, [12] resolved several contradictions and errors in the earlier papers 
mentioned above, and investigated the dynamics of the ODEs away from the 
codimension-2 point; a full description of the dynamics of the ODEs becomes 
very involved. Here we extend this work and show that small deviations from 
the exact 1 : 2 situation result in additional long- wavelength instabilities of 
otherwise stable spatially periodic patterns ('spatial quasiperiodicity'). 

The most unusual part of the dynamics of the 1 : 2 resonance problem is the 
occurrence of a structurally stable heteroclinic cycle. The local information 
that we are able compute near the equilibria on the cycle enables analysis of 
the stability of the cycle to long-wavelength spatial disturbances. 

The equations that we derive to capture these modulational instabilities are 
nonlinear PDEs, and as such, a full analysis is a daunting task. Indeed, even 
a complete classification of possible solutions is far beyond the scope of this 
paper. This paper considers those parts of the problem that can be treated an- 
alytically rather than presenting a superabundance of numerical results. Our 
analysis explores the stability of spatially-periodic equilibria and travelling 
waves, identifying those new instabilities that arc due to the inchision of mod- 
ulational effects. The dynamics near these new instabilities can be described 
by 'normal form' equations (simpler PDEs whose structure is often prescribed 
by symmetry requirements). Although the algebraic expressions may become 
lengthy, it is possible to reduce the original nonhnear PDEs to these 'nor- 
mal forms' explicitly via adiabatic elimination. Such a reduction, although 
not completely rigorously justified, can be extremely useful, particularly near 
secondary codimension-2 points (intersections of lines of codimension-1 bifur- 
cations away from the initial bifurcation point). The paper illustrates this idea 
with two explicit detailed examples (the points L and X in figure 1); in nei- 
ther case is the 'normal form' equation the Ginzburg-Landau equation. The 
dynamics near L have been well-studied in the literature, but those near X 
are more novel. 

The original study of Part 1 was motivated by a particular two-layer thermal 
convection problem, and we refer the interested reader to that paper for de- 
tailed discussion of the physical background. Wavenumber interactions in the 
ratio 1 : 2 also occur in two-dimensional thermal convection in a single layer, 
for example in the asymptotic long- wavelength equation discussed by Cox [7], 
and in non-Boussinesq convection as discussed by Mercader, Prat & Knobloch 
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[18,21]. These papers discuss situations that have one instabihty mechanism; 
by looking for solutions that are spatially-periodic with a given wavenumber 
k they locate co dimension- two points {Rc, where two periodic patterns in- 
teract, with wavenumbers in the ratio 1 : 2. Since, in these problems, Rc is 
always greater than the minimum value required to drive convection in a for- 
mally infinite layer, the results are not directly applicable to the infinite layer 
case. In contrast, this paper (formally) attempts to analyse the infinite-layer 
situation in the case where there are two distinct instability mechanisms oc- 
curring for very similar values of the control parameters. The large horizontal 
extent of the layer enables periodic patterns to be destabihsed by 'sideband' 
instabilities, for example the Eckhaus instability. We remark that a similar 
study, for 'weak' resonances, has been performed recently by Higuera, Riecke 
& Silber [11]. Their results, although very different in detail, have points of 
similarity to those presented here, for example the existence of solutions in 
the form of locahsed structures. 

The paper is organised as follows. In section 2 we derive the model equations 
and make general remarks about the dynamics. In section 3 we briefiy sum- 
marise the dynamics of the model ODEs in the absence of the modulational 
terms. Here and throughout the rest of the paper we use a single combination 
of coefficients that were used in Part 1, for ease of comparison. Sections 4 
and 5 discuss in detail the new instabilities of spatially-periodic states that 
occur when the modulational terms are added to the model. In section 6 we 
discuss the stabihty of the robust heteroclinic cycle. Analytic work shows the 
existence of a long-period periodic orbit lying close to the cycle, in quantitative 
agreement with numerical results. Section 7 briefiy highlights the coexistence, 
over a substantial region of the parameter plane, of stable non-modulated trav- 
elling waves and complex spatiotemporal behaviour. We conclude in section 8. 



2 Model equations neeir 1 : 2 resonance 

Consider a horizontal two-dimensional layer, or layers, of fiuid in the domain 
(— oo, oo) X [0, 1], using co-ordinates (x, z), i.e. of finite vertical extent but ex- 
tending to infinity in the horizontal direction x. We suppose that there is an 
x-independent state which may become linearly unstable to either of two com- 
peting steady-state instabilities, with wavenumbers nearly in the ratio 1 : 2. 
By way of illustration, the analysis of Part 1 was concerned with a two-layer 
thermal convection problem where these two instabilities corresponded to the 
onset of convection predominantly in either the upper or the lower layer sep- 
arately. The most interesting dynamics can be captured by the distinguished 
limit in which the deviation of the ratio of the critical wavenumbers from the 
exact value | is of the order of the square root of the deviation of the bifur- 
cation parameter (in this case the Rayleigh number R) from its critical value 
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Re- If wc set R — Rc = £^ then we may write the critical wavenumbers (those 
associated with local minima in the value of the critical Rayleigh number) 
as k + sq and 2k. A suitable ansatz for small- amplitude solutions near the 
codimension-2 point where the conditions for instability coincide is then 



u{x,z,t)^e [A{X,T)fi{z)e 



ia;(l+£g) 



+ B{X,T)f2{z)e'''' + c.c. 



where lengths have been rescaled so that k — 1; X — ex and T — sH are 
long length and time scales, c.c. denotes complex conjugate, and the functions 
/i 2(2;) give the vertical structure of the eigenfunction corresponding to each 
mode of instability. The amplitudes A(X, T), B{X,T) are complex-valued. 
The initial homogeneous state is symmetric under the Euclidean group E{1) 
generated by the reflection and the translations ts : x ^ x + S. 

These symmetries induce the following transformations on the amplitudes A 
and B : 



m^: x^-x: {A,B)^{A,B), (2) 

Tg: x^x + 5: {A,B) ^ {Ae'^^+"'^^,Be^'^). (3) 

The combination A'^Be^^'^^ is found to be the lowest-order translation- invariant 
combination that is not a product of the usual terms j/lp and 

The resulting amplitude equations (ignoring terms of order higher than three 
in A, B and dx) take the form 



i = A[Ati - ail^l" - + agiSe"'^^ + a^Axx, (4) 

B ^B[n2 - a2\B\^ - b2\A\^] + MV^'"^ + cBxx, (5) 

where the coefficients aj, bj and c are forced to be real by the reflection sym- 
metry (2), the dots denote derivatives with respect to the slow time scale T 
and Hi, H2 are bifurcation parameters. We remark that (4) - (5) contain both 
quadratic and cubic terms in the amplitudes A and B. To ensure a rational 
scheme of approximation in the limit e ^ we should arrange that a^, 63 
scale as e. This can be achieved in several ways, for example where a further 
symmetry which would leave the equations invariant under the sign change 
{A, B) {—A, —B) is weakly broken. 

The form of equations (4) - (5) shows that when |g| becomes large the spatially- 
averaged contribution from the quadratic terms decreases to zero. In this 
limit we formally recover the usual 'Landau' equations describing two coupled 
modes of instability in the absence of spatial resonance. This is in agreement 
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with our ansatz (1); the system is far from the 1 : 2 mode interaction point 
when \q\ ~ 0{l/e). 

Subsequent calculations are made considerably easier if we make the change 
of variable A = Ae'^*^^ to remove the exponential factors, and rescale the 
variables A, B, T and X to set — — 1 and 63 = ±1. Dropping the carat 
on A we obtain 



A = A[ni-q^ -a,\A\^ -b,\B\^]+AB + 2iqAx + Axx, (6) 
B = B[ii2 - asl^l' - 621^1'] ±A^ + cBxx, (7) 

which is the form of the equations that we will use in what follows. These 
equations have the same structure as the ODEs derived in Part 1 with the 
addition of terms giving modulation on the long lengthscale X. The term 
2iqAx captures the effect of the departure from exact 1 : 2 resonance. 

Writing A = R{X,T)e'^^^^'^^ and B = S{X,Ty'^^^^'^\ the evolution equa- 
tions (6) - (7) become 



= -q^- aiR^ - hS^] + RS cosX - ^qROx + Rxx - i?(^x)^(8) 
Re^RSsmX + 2qRx + 2Rxex + R0xx, (9) 
S^S\pi2 - a2S^ - h2R^] ± i?'cosX + c[Sxx - (10) 
S4> = ^R^ sin X + c[2Sx(t>x + S(t>xx\ , (11) 

where X = (j) — 29. When the modulational terms are omitted, we can express 
the dynamics in terms of only the two moduli and the one phase difference X 
(i.e. a reduction to a third-order system), as was done in Part 1. This leads to 
the ODEs 



R = R[fii - q^ - aiR^ - hS^] + RS cos X, (12) 
S^S[i22-a2S^ -b2R^]±R^cosX, (13) 

X^ir^-^S^sinX. (14) 

However, in the presence of modulational terms, each of the individual phase 
variables 6 and is dynamically independent. Vyshkind & Rabinovich [24] 
introduced the new variables u = S cosX, v = S'sinX (also used by Porter 
& Knobloch [20]) to produce a third-order system avoiding the co-ordinate 
singularity, but this change of variables offers no simplification in the PDE 
problem. 

The choice of the sign of the ±/l^ term in (7) has a huge effect on the dynamics. 
It was shown in Part 1 that in the case for the non- modulated problem 
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the dynamics are much less interesting than those may occur in the '— ' case. 
Indeed in the '+' case there is a Lyapounov functional V{R,9, S,4>) for the 
dynamics when 2bi = 62: 

V^(^^R^ + ^S' + hR'S' - ^S' - ii.R'' + R\ + R'{9x + qf (15) 
-R^ Sco&X+^[Sl + S^4)\]J , 

where (■ ■ ■) denotes a horizontal average. Then a direct calculation shows that 
y = - (2R^ + 2R^e'^ + ^2 + 5^0^) < 0, 

i.e. the system evolves monotonically towards a steady state. More complicated 
dynamics, for example temporal oscillations, are not possible. It follows that, 
at least when 2hi and 62 are not widely different, or the amplitudes arc small 
( equivalent ly, near to the codimension-2 point), we expect solution trajectories 
to tend asymptotically to equilibria after long times. 

For the remainder of the paper we consider the '— ' case, choosing the minus 
sign in (7). In this case it is not possible to construct a Lyapounov functional, 
and even in the absence of modulational terms the related ODE problem (12) 
- (14) can indeed display oscillatory dynamics; moreover there is a region of 
parameter space where a robust heteroclinic cycle exists and is stable. This 
cycle was analysed in Part 1 and by Armbruster et al. [1]. The parameter c 
in (7) corresponds to differential diffusion rates of the two amplitudes. Vary- 
ing c away from unity leads to Turing instabilities, well-known in the context 
of reaction-diffusion equations. We remove this complicating consideration by 
setting c = 1 in nearly all of what follows. The occurrence of Turing instabih- 
ties in mode interactions is common to all cases of strong spatial resonance and 
is of a different type to the instabilities we discuss here, since it occurs in the 
spatially-extended but exactly-resonance case. A full discussion of this second 
class of instabilities is given in a companion paper ([22], in preparation). 



3 Non-modulational dynamics near = H2 — 

In this section we summarise the relevant parts of the bifurcation sequences 
observed in Part 1 near /xi = /i2 = in the analysis of the ODEs (12) - (14), 
setting q — 0. As remarked on earlier, even in the absence of the modulational 
terms, a full analysis of (6) - (7) is extremely complicated. 

The trivial equilibrium A = B = is stable in the quadrant /ii < 0, 112 < 0. 
Near the codimension-2 point at iJ-i — 1^2 — 0, the ODEs (12) - (14) support 
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simple non-trivial equilibrium solutions of three types. The first type is a pure 
mode solution P, of the form A = 0, \B\'^ = fi2/0'2- The continuous symmetry 
of the problem implies the existence of a group orbit of equilibria; that is, the 
phase of B is arbitrary due to the underlying translational symmetry. Within 
the subspace Fix(ma;) where A,Be'R there are two equilibria, denoted P±. 
The pure mode solutions bifurcate from the trivial solution and exist when 
/i2 > 0. They are stable for fii sufficiently negative. The other equilibria are 
mixed mode equilibria of two kinds, M±, corresponding to X taking the values 
X = and X — tv, respectively. The M± amplitudes R — Ro and S — So are 
given by solutions of the equations 



where the sign choices select either M+ or M_. In fact there is a group orbit 
of each of M+ and M_ equilibria also, since although X is fixed at or tt 
respectively, there is a free choice of one of the underlying phases 9 or (f). 
When /xi is increased, holding //2 > fixed, stable M+ equilibria are created 
in a bifurcation at which the pure mode solutions lose stability. 

The M+ solution then loses stabihty (as Hi is increased further) either through 
a symmetry-breaking drift bifurcation which produces travelling waves (TW) , 
or through a Hopf bifurcation to standing waves (SW). The TW solution 
resembles the mixed-mode M_|_ in form, but it drifts along the group orbit of 
mixed-mode solutions as time evolves. The TW bifurcation is clearly a phase 
instability rather than an amplitude instability. Because the phase variables 
9 and evolve at constant rates such that X is constant and sinX 7^ 0, the 
TW solutions appear as equilibria in the {R, S, X) variables. In other words, 
the use of the {R, S, X) variables identifies all points on the group orbit of M4. 
solutions as a single equilibrium, and, in these co-ordinates, information about 
drift around group orbits is lost. 

In contrast, the SW bifurcation is an amplitude-driven instability and pro- 
duces periodic orbits that lie within the subspace X = 0. Typical curves along 
which M+ undergoes bifurcations to TW or SW solutions are shown in fig- 
ure 1. For a large region of parameter space these two curves intersect at a 
codimension-2 point, labelled A in figure 1. Because the eigenvectors corre- 
sponding to the TW and SW bifurcations are orthogonal the codimension-two 
bifurcation corresponding to simultaneous instability is a pitchfork-Hopf bi- 
furcation (in the {R,S,X) co-ordinates). 

As we increase fii for small positive fi2, the SW instability occurs first. Within 
the subspace Fix{mx) the periodic orbit created in the SW bifurcation grows 
until it collides simultaneously with the origin and the P+ equilibrium. After 
this global bifurcation the periodic orbit disappears, but now the unstable 



0^fXi±So-aiRl-hSl 

= //2'S'o ^ Rq — a2<S'o — 62-Ro'S'o- 



(16) 
(17) 
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Fig. 1. (a) Instability boundaries in the {ni,iJ,2) plane for the illustrative coefficient 
choices 0,1 = 1, 02 = 5, 6i = 2, 62 = 0, c = 1, g = 0.2. Unstable and stable solutions 
in each region are indicated with and without parentheses. The solid curve contain- 
ing the point L denotes the boundary of existence of M+. The solid curves TW and 
SW denote instabilities of M+ to travelling and standing waves respectively. The 
codimension-2 points A, X, Y, L, S and T and the various dashed and dash-dotted 
curves are discussed in the text, (b) is an enlargement of (a). 
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manifold of P+ tends asymptotically to the pure mode solution P_. By sym- 
metry, within the invariant subspace Fix{mx o r^r) the unstable manifold of 
P_ tends to P+, and a heteroclinic cycle is formed. The hctcroclinic cycle is 
structurally stable due to the existence of invariant subspaces, within which 
each connecting trajectory lies. Hence it exists for an open set of values of the 
coefficients and bifurcation parameters. Moreover, this cycle is attracting for 
an open interval of values of /ii. 

At larger /ii this cycle ceases to attract nearby trajectories as it undergoes a 
resonant bifurcation, resulting in a global loss of stability, and creating mod- 
ulated waves (MW). MW are destroyed in a Hopf bifurcation from the TW 
which themselves cease to exist in a bifurcation with the M_ states. More- 
over, several other heteroclinic cycles are possible for other combinations of 
coefficients, and at larger positive /i2- These have been investigated in detail 
by Porter & Knobloch [20]. 



4 Modulational instability of the pure mode 



Having summarised the behaviour of the system in the absence of spatial mod- 
ulations, we now allow q to take non-zero values, and examine the possibility of 
new modes of instability due to the spatial frequency mismatch. We first note 
that the results of Part 1, summarised in the previous section, apply within the 
subspace of non-modulated (X-independent) solutions. To examine the sta- 
bility of the pure mode solution P_|., given by A — and B — Bq = 112/ o.'i-, 
we substitute the ansatz 

A{X, T) = ai(r)e^^^ + a2(T)e-^^^, 

B{X, T) = So(l + /3i(r)e^"^^ + P2{T)e-'^''), 

into (6) - (7) and linearise. The linearised dynamics for ai^2 and /3i,2 decouple 
and P+ is found to be stable to the perturbations /3i,2 for all wavenumbers m. 
The linearised system for the perturbations ai 2 is 



fii-q^- hiBl - 2qi - f Bq 

Bo iii-q^ -biB^ + 2qi-i\ 

~ (18) 

This matrix has trace tr(^^) and determinant D{P): 
tT{e')=2{^,-q'-f-hB',), 



1 






[02) 
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Die'') = {^,^ -q^-f- hBlf - AqH'' - 



No oscillatory bifurcation is possible since tr(£^) = implies D{i'') < 0. 
A steady-state instability to perturbations with wavenumber i occurs when 

D{f) = D'{f) = since D"{f) = 2 > 0. These conditions show that the first 
instability of P+ may be to perturbations either with £ = or with i non-zero. 
When H2 > 0'2{^^i + Q^)/bi the first instability is to £ = and occurs along 



the curve jdi = 61/^2/02 + q'^ — YA*2/a2- When /X2 < a2(/xi + q^)/hi the first 

instability is to finite £ = £c > and occurs along the line 



These two instability curves meet at the point (/Xi,/^!) — {'^^iQ^ ~ ) '^0'2<i^) 
marked L on figure 1. The gradients of these curves are equal here, so the 
transition between instabilities is smooth. Along the line OL (from the origin 
to L) the most unstable wavenumber, l^i is given by £^ = /ii — 61/X2/Q2 + (fi 
i.e. < < and Ic increases monotonically as the origin is approached. 

At any point on the interior of the line OL we can fully describe this bifurcation 
by the usual Ginzburg-Landau equation, since for the PDEs (6) - (7) this is an 
instability of a uniform state to a nonzero-wavelength perturbation, and the 
growth rate of a perturbation with a wavenumber far from £c is negative and 
bounded away from zero. We have carried out a weakly nonlinear perturbation 
expansion near OL to investigate whether this bifurcation is subcritical or 
supercritical, details of which are given in Appendix 1. The resulting analytic 
expression is cumbersome, but it is possible to deduce that the bifurcation 
is always supercritical when ai > is large enough. This calculation can be 
carried out analytically without assuming c = 1; if c is large compared to 
unity and ai is small it is possible for the bifurcation to be subcritical. For the 
illustrative set of coefficients used in figure 1 the bifurcation is supercritical 
along the whole of OL. 

The dynamics in a neighbourhood of the codimension-2 point L cannot, though, 
be described by the Ginzburg-Landau equation, since the instability wave- 
length £c tends to zero as L is approached. A codimension-2 point identical in 
structure to L occurs in the analysis by Coullet & Repaux [5] of a pattern- 
forming instability subjected to an external nearly-resonant periodic forcing. 
They term L a 'Lifschitz point'. Through asymptotic expansions near L it is 
possible to describe the dynamics in terms of a single 'normal form' equation 



which governs the evolution of a small perturbation A{^, r) e R to the P+ 





dA 



(19) 
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solution; ui and U2 are new bifurcation parameters, defined so that the point 

L corresponds to z/i = z/2 = 0, r is a new scaled time variable and ^ is a new 
long spatial scale associated with the smallness of £c- Equation (19) is known 
as the Extended Fisher-Kolmogorov equation and it's properties have been 
extensively investigated ([3,19,2]). It is easily shown that there is a Lyapounov 
functional for the dynamics, and so there arc only steady state solutions at 
long times. However, these states need not be periodic in space; and in fact the 
solutions can have very complex spatial structure when i/i > 0, corresponding 
to the region ni > hiB^ — Bq + q^. 



5 Modulational instabilities of the mixed modes 



M_(. undergoes a plethora of different bifurcations. In this section we will dis- 
cuss three new instabilities to long-wavelength disturbances that occur. We 
examine a co dimension- two bifurcation where two of these curves meet. We 
also discuss other codimension-two points that occur where one of these long- 
wavelength instabilities meets a bifurcation curve from the non-modulated 
problem. Our discussion is organised by the sequence in which these bifurca- 
tions appear in figure 1 as //i increases. 

Let A = Ro{l + aie'^^ + ase"^^^) and B = So{l + P^e'^^ + /^ze"'^^) where Rq 
and y^o satisfy 



= fii-q^ + So-a,Rl-hSl (20) 
= fi2S-Rl- a2S^ - b2RlSo. (21) 

After substituting into (6) - (7), linearising and changing to the sum and 
difference variables Q!± = cti ± 0:2, /?± = /?! ± /?2, we obtain the linearisation 
matrix 



^ Roa+ ^ 


/ 


-2aiRl - 


Ro( 


1 - 26150) 


-2g£ 


\ 










-2Ro{l + b2So) 




2a2S^ - ce 










s,p+ 


RqCX- 




-2qe 







-25'o - £2 


Rq 




RqU^ 




V 










-2Ro 


So " / 




\ SoP- J 



(22) 



The characteristic polynomial of this matrix can be written as -P(A) = A"^ + 
i(£2)A3 + 5(£2)A2 + C{e^)\ + Z)(£2). Note that P(A) always has a root A = 
when i = 0, hence we may write Z)(£^) = i'^E{i'^) where Ei^i"^) is a cubic 
polynomial. This is due to the underlying translation symmetry. Steady-state 
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instabilities at £ = > occur when E{ll) = E'{il) = and E"{il) > 
0. Oscillatory instabilities occur when A = a{P) + ic<j(£^) satisfies a{il) = 
a'{il) = 0, a"{il) > and iu\il) > 0, where a{e) and uj{e^) are real. These 
requirements yield the conditions 

C"^ - ABC + A'^D = 0, 
{AC' - CA'){2C - AB) - A\B'C - Ab') = 0, 

C/A > 0, 

for an oscillatory bifurcation with > 0. Both the steady-state and the 
oscillatory cases give two conditions; in conjunction with (20) - (21), these 
conditions enable the determination of bifurcation lines in the (//i,/X2) plane. 
Clearly, a steady-state instability at = occurs when £"(0) =0 and £"(0) > 
0; similarly an oscillatory instability occurs when a{0) = and a'(0) < as 
long as a;^(0) > 0. In terms of the coefficients of -P(A) these conditions are 

(7^ - Abc + A^b = 0, 

(Bc -Ab- cyAy 
- — - — — — ^ — ^ <o, 

2{4:D-AC-B^) 

C/A>0. 

5. 1 Instability of M+ near the Lifschitz point L 

The first new instability of M_|_ that wc find is a steady-state bifurcation to 
spatially-modulated solutions, i.e. the most unstable wavenumber is non-zero. 
Tills occurs along the dashed curve LT emanating from the Lifschitz point 
L. This instability is part of the generic bifurcation structure near a Lifschitz 
point and hence the existence of this curve can bc shown by analysis of the 
normal form (19). We have followed the bifurcation curve to larger values of 
//i where it becomes asymptotic to the line //i = g'^ as //2 — oo. Near L, when 
Oi is large we expect the bifurcation to be subcritical by comparison with the 
results of Coullet & Rcpaux [5]. This bifurcation curve meets a second such 
non-zero-£ instability curve (this one starting from the point §) at the point 
T on figure 1. At T there is a Hopf/steady-state mode interaction coupled 
to a phase mode which has zero growth rate at zero wavenumber; figure 2 
shows the variation of the real parts of the eigenvalues A(£^) at this point. 
The dynamics near the codimension-two point T depend strongly on the ratio 
of the wavcnumbers involved in the instabilities. This ratio varies considerably 
with g, as shown in figure 3. For the illustrative coefficient set used in figure 1, 
when q — 0.097 there is a codimension-three bifurcation as § coUides with T. 
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0.005 




Fig. 2. Codimension-two instability of at the point T, 

{l_ii,fj,2) = (0.00929,0.0149) for q = 0.2 and the coefficients of figure 1. Real 
parts of the eigenvalues A of the linearisation about are plotted against 

wavenumber £'^. Solid lines denote real eigenvalues, dashed lines give the real part 
of complex conjugate pairs. Note the mode with zero growth rate at ^ = due to 
the translational symmetry of M+. 
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q q 

Fig. 3. (a) Critical wavenumbers for the steady-state (solid line) and oscillatory 

(dashed line) instabilities at T (see figure 1) as a function of the wavenumber mis- 
match q, in the range 0.097 < q < 0.55. (b) Ratio of the critical wavenumbers in 
(a) as a function of q. 

As we approach this collision, the wavenumber of the oscillatory instability at 
T goes to zero and so the wavenumber ratio is formally infinite there. 

The critical wavenumbers for the steady-state instability can be seen from 
figure 3(a) to be slightly greater than the mismatch parameter q, in contrast 
with the modulational instability of the pure mode solution P+ discussed in 
section 4 where the maximum critical wavenumber of instability is exactly q. 
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It turns out that for M+ the maximum instabihty wavenumber is 2q. This can 
be easily shown by examining near the hne = . fj.2 > 0. When //i — 
is small and negative the amplitudes Rq and 5*0 for M_|_ are, solving (20) - (21) 
to leading order: 

ai/i2 — 1 

/ii - g2 

•^o— r- 

aiA*2 - 1 

Then, substituting these expressions into the term E{^P) from the character- 
istic polynomial P{\) we obtain 

E{e) = - - V) + o{\t,, -q'\). 

The double root at £^ = ^2 comes from the LW instability; it is the root at 
^2 — ^-^^^ comes from the continuation of the steady-state instability curve 
through T towards /xi — ~ 0. From numerical investigations we conjecture 
that the wavenumber of the instability evolves monotonically along the curve, 
but there does not seem to be a straightforward way to verify this analytically 

5.2 Instability of to standing waves 

The dash-dotted hne TS in figure 1(b) marks the instabihty boundary of 
M_|_ to the second new bifurcation involving spatial modulation; an oscilla- 
tory instability to temporal oscillations (standing waves) with non-zero spatial 
wavenumber. The dynamics near the point § is analogous (but time-periodic 
rather than steady) to that near the Lifschitz point L, and can thus be de- 
scribed by a similar extension of the complex Ginzburg-Landau equation: 

dA 

— = (i/i + iuoiui, U2))A - (1 + i(3)A\A\^ + {u2 + ia)% - (1 + 

(23) 

governing the evolution of a complex- valued perturbation A{^, r) to the 
solution. As before, ui and 1/2 are new bifurcation parameters; the point § 
corresponds to i/i — 1/2 — 0, cuq is the non-zero frequency of the instability, a, 
f3 and 7 are real parameters, r is a new scaled time variable and ^ is a new 
long spatial scale. Clearly the dynamics of (23) are at least as complicated as 
those of (19). The wavenumber of this oscillatory instability increases as ^2 
decreases along the curve (from S to T) . 
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5.3 The long wavelength phase instability o/M+ 



The third new instabihty occurs, at larger /i2, along the curve LW in fig- 
ure 1(a). This instability is a long- wavelength steady-state bifurcation; as q 

is increased from zero, the LW curve splits off from the TW curve along its 
entire length. Since, in the absence of modulational terms, M_|_ are unstable 
first to SW below A and unstable first to TW above A, the LW instability (at 
least for this combination of coefficients) is the instabihty of M+ that occurs 
first in an intermediate range of /i2, between the new codimension-two points 
X and Y. The LW and SW curves intersect at Y and the LW and TW curves 
intersect at X, see figure 1(a). Since the LW instability curve coalesces with 
the TW instability as g — > it is clear that the LW instability must also be a 
phase instability rather than an amplitude instability. 

It is important to observe that the point A, where the TW and SW curves 
cross, is now 'shielded' by the LW curve. The dynamics near A were inves- 
tigated in detail in this context by Julien [12], and in greater generality by 
Landsberg & Knobloch [15]. In the presence of modulations the dynamics near 
A are expected to be less relevant to the observed dynamics in a spatially- 
extended physical system. However, it is worth noting that by varying the 
diffusivity ratio c away from unity it is possible to make the points X and 
A coincide; this would result in a codimension-three bifurcation involving the 
TW, LW and SW instabilities of M+. 

Reduced descriptions of the dynamics near X and Y can be derived from 
the PDEs starting either from (6) - (7) or, more conveniently here, from the 
modulus/phase representation of the dynamics given by (8) - (11). In this 
subsection we will present the reduction near X in some detail, and comment 
only briefiy on the dynamics near Y. 

Near X we strive to eliminate the equations (8) and (10) for the evolution of 
the moduli R and S by adiabatic elimination to leave a pair of real equations 
for the phases and 0. It turns out to be more convenient to describe the 
dynamics in terms o{X = (j) — 2d and 9^ since the TW instability involves only 
the combination X. From (14) it is clear that the TW instability occurs when 




(24) 



where Rq and 5*0 satisfy (20) - (21). The LW instability occurs when £'(0) = 
where £'^E{i'^) is the determinant of the linearisation matrix (22): 



^(0) 




'0 



{2cS^o - ^o)(2(aia2 - hb2)Sl + (62 - 2b,)S^o 
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-a^Rl + So) + 2q\2a2Sl - RD] . (25) 

Note that when c = 1 and g = this condition reduces to the condition for the 
TW instabihty (24). When c—1 the algebra simphfies substantially (and so 
we set c = 1 for the remainder of this section); solving (24) and (25) together 
yields the moduli Rq = \/2/a2, Sq = 1/ a2 the point X where the LW and 
TW instabilities coincide. Substitution into (20) - (21) gives the co-ordinates 
of the point X in the (/xi, 112) plane: 



(/^f, l4) = ((2ai + 61 - a2)/al + q\ (802 + 262)/a^). 

The adiabatic elimination of R and S at the point X proceeds in the usual 
manner, but the scaling for X differs from that which might be expected. Our 
choices of scalings are determined completely by the requirement to balance 
the hnear terms in the reduced equations. This balance then introduces two 
specific nonlinear coupling terms at leading order. We write 

R^Ro + sVo + sVi H , X^eX, (26) 

S^So + e^so + eh-, + • • • , T = s^T, (27) 
{eA)=e{l~4>), X^s^X, (28) 

//I + II2^ l4 + ^^f^2. (29) 

Substituting into (8) and (10) and dropping the tildes gives, at 0(£^): 

2aiRlro + {2biSo - l)RoSo = Rof^i - IqRodx, 
2i?o(l + &2'S'o)ro = Sojji2- 

Hence 



'^°=2fi„(l + W„)' 

2qex ^0^2 .01 N 

'° 2hSo-l {2hSo-l){l + b2So)' ^ ^ 

Substituting the scalings (26) - (29) into the 9 equation (9), dropping the 
tildes and cancelling a factor of yields 



Roe^ R^SoX + 2qrox + RoOxx + O {e^) . (32) 
Prom (30) we see that rox = and so (32) simplifies to give 

e^SoX + 9xx + 0{s''), (33) 
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which is the first of our pair of reduced equations. It turns out that we need to 
compute the term ri to determine the leading order evolution of X correctly. 
From (10) at 0{e^) we find 



2Ro{l + 62'S'o)ri = sofi2 - 3025*050 - (1 + hSo^Q - 262-Ro?^o-so 

+SOXX - SQ{(t)xf - So, 

which, after substituting (30) and (31), becomes 

12ga2-S'o[(l + b2So)fii - SofJi2\ 2g/i2 



2Ro(l + b2So)ri 



(2hSo - 1)2(1 + b2So) (1 + b2So)(2b^So - 1), 



ASo{9xf + ^, , {Ox - Oxxx) + const, 



where const denotes terms independent of X , and the last term is equal to 
2qSoXx / {2biSo — 1) at leading order, using (33). 

Now we turn to the (unsealed) X equation, formed by combining (9) and (11): 



X^(^^-2S^X + Xxx + ^Sx(l)x - ^Rx - ^RxOx. 

After substituting the scalings (26) - (29) and dropping the tildes, the terms 
at 0{e^) are found to be 



X- 



2Ro ( Rl\ 

Vo - 2 + So 



Q u 1 ' r<2 
'-'0 \ >-'0 



4 Aq 
X + Xxx + -^soxOx - s-nx, (35) 



using the fact that tqx = 0. After substituting for sqx and rix using (31) 
and (34) (noting that the terms in (34) indicated by 'consf do not appear) 
we obtain 



X^CoX + (1 + ^4)Xxx + iiOxx + ^2X9x + ^sOxOxx + 0{s^), (36) 
where 



6 = 



At2 



0.2 



4a9 



02 + 62 (26i - a2)(a2 + 62) 



402/^1 
2bi — a2 ' 



-12q^aljli 



(a2 + 62)(26i-a2)2 
8ga2 



+ 



2q^alfi2{Q + 2bi - 02) 
(02 + 62)2(261-02)2 ' 



26i — a2 ' 
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_ 4qal{2h-3a2-2b2) 
(a2 + 62)(26i-a2) ' 
-2q^al 
(a2 + 62)(26i-a2)' 

One simple consistency check is that coefficients containing odd powers of 
q multiply terms with odd numbers of X-derivativcs (and likewise for even 
powers of q). For the illustrative coefficient choices ai = 1, 02 = 5, 61 = 2 and 
62 = we obtain the leading order reduced equations 



X = (20/ii - 3/i2)X + (1 + 50q'^)Xxx 

-1250g2(6/ii - il2)exx - ^OqXex + 220g^xfe, (37) 

0^lx + 9xx. (38) 
o 

Non-modulated M+ states correspond to X = 0, 6* = 6*0 as there is a circle of 
equivalent M+ states related to each other by spatial translations. We consider 
the state ^0 = without loss of generahty. We now use (33) and (36) to 
examine the two distinct linear instabilities of the state X = 6 = that are 
possible. Substituting X = X{T)e'^^ and 9 = 9{T)e'^^ into (33) and (36) and 
linearising we obtain the Jacobian matrix 



Co -^'(1+^4) -^1^ 

1/02 -e 



J- 



which has trace and determinant 



detjif) + U) + ^'(6/02 - ^0). 

Hence detj'(O) = ^i/a2 — Co- The bifurcation to TW occurs when = 0; it is 
the initial instabihty of M+ when ^1/02 > 0, i.e. (for the illustrative coefficient 
set) when 20/ii — 3/i2 = and 6/ii — jj,2 > 0- Similarly, the LW instability 
occurs when ^1/02 = ^0; it occurs before the TW instability if ,^0 < 0. These 
conditions become (20 + 1500g^)/ii = (3 + 250g^)/i2 and 20/ii < 3fi2 for our 
coefficient set. These lines are sketched in figure 4 and clearly correspond to 
the behaviour of the TW and LW curves near X in figure 1(a). 

The form of the equations (37) - (38) stems directly from the requirement 
of invariance under x-reflection (2), which sends {X,9,dx) {—X, —9, —dx)- 
We note that the TW bifurcation is one of the secondary instabilities of spa- 
tially periodic patterns classified on symmetry grounds by CouUet & looss 
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Fig. 4. Sketch bifurcation lines near the codimension-two point X where the TW 
and LW curves meet, for the illustrative coefficient choices ai = 1, a2 = 5, 61 = 2, 
62 = 0, c = 1, g = 0.2. M+ are stable above and to the left of the LW and TW 
curves. The TW bifurcation (solid/dotted line) occurs along the line 20/ti — 3/x2 = 0; 
it is the initial instability of Af+ when 6/ii — /i2 > (indicated by the solid line). 
The dotted part of the line indicates instability in the absence of long-wavelength 
modulations. The LW bifurcation (solid/dash-dot-dot-dotted line) occurs similarly 
when (20 + 1500q^)fli = (3 + 250q^)fi2 and 20/ii < S/ia. The solid part of the LW 
line indicates that part of it for which it is the initial instability of M+. 

[4]; their equations (8a,b) closely resemble (37) - (38), although theirs con- 
tain only one bifurcation parameter since they are concerned with classifying 
codimension-one instabilities. By including a second bifurcation parameter, 
K,2, wc are able to capture the codimension-two transition between the TW 
and LW instabilities. 



5.3.1 Eckhaus instability dynamics near the LW bifurcation 

At the codimension-one LW bifurcation, instability to modes of arbitrarily long 
wavelength occurs; near this bifurcation we may adiabatically eliminate X to 
derive a single real equation describing the dynamics. The relevant scalings 
are 

e = ee, X = eX, 
X = e% f = e% 

we introduce a bifurcation parameter u by writing e'^u = 1 — Ci/('22Co)- On 
substituting these scalings into (36) we obtain 
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CoX — —CiOxx + £^ 



-yxxxx + — : ^3 yx^xx 



6 V Co 

and we substitute this expression for X into the 9 equation (33) to obtain 



ti — i/t>xx -\ T2 — i^xxxx -\ T2 t7xx , i -jy j 

02^0 Ct2^o 



to leading order. Equation (39) is identical in form to that which describes 
the weakly nonlinear behaviour of the Eckhaus instability. Like the Eckhaus 
instability, the LW bifurcation is therefore always subcritical. 



5.3.2 Instability of travelling waves near the TW bifurcation 

In a similar way, the dynamics of travelling waves, near the TW bifurcation, 
near X, can be investigated. At the TW bifurcation, 9 can be eliminated 
(again, adiabatically) from (37) - (38) and the resulting single real Ginzburg- 
Landau equation for X describes the bifurcation leading to TW solutions. The 
scahngs leading to (33) and (36) cannot be chosen to include a term X^ in (36) 
which would be required to capture stable finite-amplitude TW states near 
the bifurcation point. In the absence of modulations a different set of scalings 
can be chosen which arc able to include this term. In this way the sub- or 
supercriticality of the TW bifurcation can be easily computed. 

However, it transpires that spatially-periodic TW arc unstable to modulational 
disturbances and this instability is indeed captured by the reduced equations 
near X. To illustrate this, we compute the stability of the X-independent state 
X = Xq constant, 9 = XoT/a2 for (33) and (36). 

Let 

X = Xo(l + ae^'"" + C.C.), ^ = ^(T + /3e^« + c.c), 

02 

then, on substituting into (33) and (36) and linearising we obtain 

+ «[6 - ^'(1 + 6)] + -(i^6Xo - ^'6), 

02 

Eliminating a yields a single linear, constant coefficient ODE for /3: 

'$ + /3(£2(l + ^4) - 6) + /9(^'(l - 6 + 6/«2 +^'(1 + 6)) - i^6Xo/a2) =6- 
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Solving the homogeneous equation for the complementary functions (3 — & 
we find 



6 - ^^(1 + ^4) ± ^J{Co - + ^4))^ - 4(£^(1 - Co + a/«2 + ^^(1 + ^4)) - i^6Xo/a2)) 



For sufficiently small ^o, the growth rate A+ is more positive than the growth 
rate of the state X = Xq, ^ — XoT'/a2 towards a fuUy-nonhnear TW equihb- 
rium, hence stable TW are not anticipated to appear close to the codimension- 
two point X. 

Near X, numerical results (see figures 5 and 6) show M+ solutions losing 
stabihty first to a TW perturbation which generates a TW-hke 'transient', 
and then the occurrence of an instability to spatial modulations. The final 
state is, however, steady, and consists of an extremely long-wavelength, but 
0(1) in amplitude, modulation - see figure 6. Throughout most of the domain 
Q and increase linearly with X, while X remains close to zero. 

From the form of the equations (8) - (11) it is clear that families of solutions 
with dx and 0x constant and non-zero are possible, corresponding to exactly 
spatially-periodic states with wavenumber close to unity. Some analytic in- 
vestigation of these solutions should be possible. It may well be possible also 
to look analytically at other classes of solution, for example homoclinic orbits 
corresponding to spatially localised structures in infinite domains, or solutions 
where R and S are spatially periodic, but the possibilities are too numerous 
to discuss further here. 

5. J). Dynamics near the codimension-two point Y 

Below the point A on figure 1 the LW and SW curves intersect at yet another 
codimension-two point, labelled Y on figure 1(a). A similar analysis to that 
near X could be carried out here, leading to a pair of equations for the SW 
instability (which leads to time-periodic variations in solution amplitude) , and 
the LW (phase) instability. Although we have not computed the reduced 'nor- 
mal form' in this case, its structure is simple to derive and we include it for 
completeness. The relevant reduced equations describing the dynamics of this 
bifurcation are coupled ODEs for two variables {z, 9) where z E C gives the 
perturbation in the direction of the SW instabihty and 9 E M. describes the 
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Expanding this expression for A+ up to 0{£^) we obtain 
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Fig. 5. Time evolution of perturbation from the M+ equilibrium near X. /xi = 0.018, 
H2 = 0.7; all other coefficients are as in figure l.The four panels show spatial averages 
of: (a) \A\ (solid line) and \B\ (dashed line), (b) \Ax\ (solid) \Bx\ (dashed), (c) X 
(solid), 9 (dashed), (d) \Xx\ (solid), 1^x1 (dashed). An initial transient growth of a 
TW mode (illustrated by the constant rate of evolution of 6 for 900 < t < 1500) 
gives way to a spatially-modulated steady state at t = 1600. 
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Fig. 6. Steady spatially-modulated state reached at the end of the numerical inte- 
gration illustrated in figure 5. /ii = 0.018, H2 = 0.7; all other coefficients are as in 
figure 1. (a) 1^1 (soUd), \B\ (dashed), (b) X (solid), 9 (dashed). 

LW instability. On symmetry grounds, the equations will be invariant under 
the transformation {z, 9, dx) — {z, —9, —dx)- We also make use of the normal 
form symmetry z — > e^'^z which appears naturally in Hopf bifurcation prob- 
lems, and we use the fact that only spatial derivatives of 9 will appear because 
the value of 9 itself is dynamically unimportant. Under these constraints the 
reduced equations (including terms up to cubic order in {z,9,dx)) take the 
form 



i = (Ai + \u)z - ioz\z? + iiz9x + (40) 
9 = -\29xx - iz^xxxx + 1^4 {zxz - zxz) Csklx, (41) 

where .Co,^i,S ^ ^ ^-iid .^3, .^4, .^5 G M are undetermined coefficients, a; > 
is the frequency of oscillation at the Hopf bifurcation, and Ai, A2 are real bi- 
furcation parameters. Due to the large number of undetermined coefficients, 
space does not permit a detailed investigation of this bifurcation here. How- 
ever, since the symmetry does not permit a linear term in z in the 9 equation 
it is clear that the dynamics of (40) - (41) are not related to those of (33) 
and (36). The codimension-one bifurcation that occurs for Ai = and A2 < 
was one of the 'normal forms' identified by CouUet & looss [4] and has been 
explored numerically by Lega [16] and Daviaud et al. [10]. In particular, these 
authors identify two distinct regimes of spatio-temporal chaos depending on 



24 



the choices of the couphng coefficients in (40) - (41). 



5.5 Instability of M_ to long-wavelength perturbations 

Finally in this section, we briefly discuss the dynamics in the quadrant ni > 0, 
//2 < 0. The spatially-periodic equilibrium state M_ exists in the whole of this 
region and is stable when is sufficiently negative, for a fixed /ii > 0. The 
amplitudes R — Rq and S — So satisfy 



and X — n. For g = 0, M_ states lose stabihty to TW solutions as //2 is in- 
creased at fixed positive /ii. These TW states have X ~ tt near the bifurcation 
since X = TT for M_. As for this phase instability also generates a distinct 
long-wavelength instability when q ^ 0. 

Following the usual linearisation of (8) - (11) about M_ we find that a long- 
wavelength instability occurs when 



{2cS^ - Rl){2S'^{b,b2 - aias) + (62 - 2bi)S^ - Sq - a^Rl) = 2q\Rl + 2a2Sl). 

(44) 



Since the bifurcation to TW occurs when R^ — 2Sq as before, it is clear that 
in the case c — 1 these bifurcation curves do not intersect and no further 
codimension-two bifurcation points appear. In the limit where we keep ~ 
0(1) fixed and let fi2 —00 we see from (42) - (43) that Sq ^ and 
Rq ~ 0(1). Then (44) is satisfied when aiR^ ~ 2g^ to leading order. From (42) 
this gives the asymptotic behaviour of the curve of LW instability: /ii ~ Sq'^, 
illustrated in figure 7. 

In summary, in this section wc have identified various codimension-one and two 
bifurcations involving modulational instabilities. In particular, two codimension- 
2 points, X and Y, organise the interaction of the new LW instability with 
the previously-studied TW and SW bifurcations. Reductions of the original 
PDEs (6) - (7) to 'normal form' equations describing the dynamics near these 
bifurcations enables us to gain insight into the nature of these bifurcations, 
despite the fact that solutions, even of the reduced equations, may have ex- 
tremely complicated spatial structure. 



= yUi - - 5*0 - aiRl - 61 5*0, 
= H2S0 + RI — 025*0 — 62-Ro'S'o) 



(42) 
(43) 
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Fig. 7. Region of stability of the mixed mode M_ for the illustrative set of coefficients 
ai = 1, 02 = 5, 6i = 2, ^2 = 0, c = 1, g = 0.2. The region of stable M_ is bounded 
by TW (solid line) and LW (dash-dotted line) instabilities. Note the asymptotic 
behaviours of the LW curve: at large negative fi2 we find jj^i ~ 3q^ = 0.12 for this 
set of coefficients; at large fii it closely follows the TW curve. 

6 The heteroclinic cycle 

One of the most interesting features of the ODE problem (12) - (14) anal- 
ysed in Part 1 is the existence of a robust heteroclinic cycle between pure 
mode solutions related by a half-wavelength spatial translation; for example 
P+ and P_. For the coefficient values selected in this paper this heteroclinic 
cycle exists for an open region of the {^1,^2) plane, abutting the origin. Its 
formation relies on the existence of pairs of two-dimensional invariant sub- 
spaces for the dynamics: for example within the two-dimensional subspace 
Fix(mj;) = {Im{A) — Im{B) — 0} the point P+ is a saddle and P_ is a 
sink and, if various other conditions are met, a connecting trajectory between 
them exists. The second connection is then forced to exist by symmetry, and 
is contained within the subspace Fix(ma; o r^-) = {Re{A) = Im{B) = 0}. 

For small positive 112 and increasing /ii the cycle is formed after a global bi- 
furcation that involves the intersection of the unstable manifold of P+ and the 
stable manifold of the origin; this bifurcation also creates or destroys the SW 
periodic orbit. At larger values of 112 the connecting trajectory appears after 
a saddle-node bifurcation marking the boundary of existence of a further two 
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M+ equilibria. Hence this curve of saddle-node bifurcations also bounds the 
region of existence of the cycle. At larger fii the cycle ceases to exist where 
the pure mode equilibria P± gain stability in a pitchfork bifurcation with M_ 
equilibria. These bifurcations proscribe the region of existence of the cycle. 
Within this region of existence, a further curve separates regions where it is 
stable or unstable. At this stability boundary a branch of modulated waves 
(MW) bifurcate from the cycle in a global ('resonant') bifurcation. Theoretical 
work by Armbruster et al. [1] and by Proctor & Jones [23] was confirmed by 
the general stability results of Krupa & Melbourne [14, section 6.1]. It turns 
out that the natural condition (that the ratio of eigenvalues in the 'contracting' 
and 'expanding' directions should be greater than one for stability) is neces- 
sary and sufficient. This yields the condition /ii < 6i/X2/a2 for the stability 
boundary of the cycle. 

When modulational terms are included, the instability of P+ to P_ that is 
necessary for existence of the cycle implies that P+ will also be unstable to 
sufficiently long-wavelength perturbations, so it might be expected that the 
cycle also could not be stable to long-wavelength perturbations. Moreover, the 
subspaces {/m(/l) = Im{B) = 0} and {Re{A) = Im{B) = 0} are no longer 
invariant for the dynamics; spatial variations of the amplitude A drive the 
evolution of the phase 9 when q 0, see equation (9). Of course, the cycle 
still exists, within the subspace of solutions with no spatial modulation at 
all. But, in large domains it cannot be asymptotically stable (that is, points 
in a full neighbourhood - in some appropriate sense - of the cycle converge 
to it) as the equilibria P± will be unstable to perturbations of a sufficiently 
long wavelength, as well as to each other within the subspaces containing the 
heteroclinic connecting orbits. 

Numerical simulations show that spatially-modulated perturbations eventu- 
ally grow, when q ^ 0, and the simulation converges to a periodic orbit rather 
than showing the characteristic increases of time spent near each equilibrium 
that would indicate convergence to the robust heteroclinic cycle. 

In this section we analyse the behaviour of trajectories close to the cycle using 
a linearised stability analysis near the equilibria Pj- and under the commonly- 
used assumption that trajectories near the cycle spend very little time passing 
between neighbourhoods of the equilibria. Our aim is to determine conditions 
for existence of periodic orbits lying close to the cycle. We restrict our attention 
to periodic orbits that spend equal amounts of time near each equilibrium since 
this feature is observed in numerical work. 

Recall from section 4 that the linear stability of P+ to perturbations ~ ai&^^ + 
A{X,T) with wavenumber i is given by the matrix (18) acting on 
the amplitudes (ai, a2). For simplicity we write this matrix as 
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where 



^ ao-q bo ^ 
bo ao + q J 



ao = //i - 
bo^Bo, 
q^2q£. 



q' - b,Bl - e 



(45) 
(46) 
(47) 



Note that a necessary condition for the existence of the cycle is that the 
eigenvalues A± of Mi, when evaluated when ^ = 0, satisfy A_ < < A+, so 
that P+ is a saddle point. The form of Mi implies that this is equivalent to 
requiring < when £ = 0, and hence in the rest of this section we assume 
Oo < 0. The corresponding linearisation around P_ is denoted M^. 



Mo 



^ao-q -bo ^ 
y -bo ao + q J 



The form of M2 is determined entirely from Mi and the equivariance of the 
ODEs (6) - (7). 

We assume that trajectories close to the heteroclinic cycle spend equal amounts 
of time T near each equilibrium, and that we may ignore the time spent 
travelling between neighbourhoods of the equilibria. Hence a perturbation 
■^0 = (tti, 02)^ near P+ evolves to the point 



vi = Svo = exp(M2T) exp(Mir)t'o, 
where the matrix S is given by 

2ble+e_ + qel{q - Co) + qel{q + Co) -boq{e+ - e_)^ 

^•^0 \ boq{e+ - C-f 2hle^e- + qe%{q + Cq) + qe^_{q - Cq) y 



and Co = \J(p + bQ and e± = exp[(ao ± Co)T]. The eigenvalues of S then 
determine the stabihty of the cycle. We are particularly interested in the de- 
pendence of these eigenvalues on the travel time T, the selected wavelength £ 
and the wavenumber mismatch parameter q. 

When q = 0, S simplifies enormously: it is diagonal with eigenvalues e+e_ — 
exp{2aoT) where ao is defined in (45). Since ao < we have contraction of 
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the perturbation vector vo under successive iterates of S. This corresponds to 
trajectories passing repeatedly through neighbourhoods of P±, and converg- 
ing to the cycle. This result is independent of the time T spent near each 
equilibrium. 

For general q ^ the eigenvalues of S will, however, depend on T. Moreover, 
for any q ^ 0, S will have eigenvalues of modulus greater than unity when T 
is taken to be sufficiently large (so that e+ is sufficiently large). So, in addition 
to the discussion above, we can conclude directly from the form of the map 
that for any non-zero q the cycle is unstable. 

Numerical simulations indicate that trajectories remain close to the cycle, 
though, so it is natural to ask whether stable dynamics close to the cycle 
is possible. One possibility is the existence of a long-period periodic orbit. 
That such an orbit might exist is motivated by the observation that S has 
an eigenvalue greater than unity for T ^ 1, i.e. trajectories very close to the 
cycle are pushed further away from it, but that for T ~ 1 the eigenvalues 
of S may lie within the unit circle, indicating that trajectories that start 
further away from the cycle move closer on successive passes near the cycle. 
This behaviour is confirmed by numerical simulations (see figure 8), where 
we choose an initial condition lying extremely close to the subspace of non- 
modulated solutions, but not close to the heteroclinic cycle. Thus the first 
few travel times T between neighbourhoods of P-t occur with T small, and 
the perturbation decays. Then, as the trajectory becomes closer to the cycle 
T increases and the perturbation grows (exponentially on average) until the 
trajectory converges to the periodic orbit. 

The above argument is unusual in that it suggests that wc can extract nonlin- 
ear information (about the period of the orbit) from a purely linear calculation 
(that which leads to the form of S). This is because we are essentially using 
T as a proxy for the closest distance between the periodic orbit and the hete- 
roclinic cycle. Hence solving for the position of the orbit and solving for T are 
really the same thing. By setting up the usual 'small box' approach we could 
construct a return map, fixed points of which would correspond to periodic 
orbits. By using T instead, we are able to circumvent the need to compute this 
return map in full in order to extract an estimate for the period. Essentially 
we construct the component of the return map in a pair of directions orthogo- 
nal to the invariant plane corresponding to unmodulated solutions. Then the 
condition for locating a fixed point of this part of the map is identical to the 
condition that S has an eigenvalue of 

At least for small q and large times T we expect that the behaviour outlined 
above could be captured by the map S based on linearisation near P±. We 
compute the approximate period P{q) ~ 2T(g) of a stable periodic orbit by im- 
posing the condition that the larger eigenvalue m+ of S is (the other (m_) 
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Fig. 8. Temporal evolution of the first few Fourier modes (writing 

A{X,T) = YljLo^j'^^'''^^'^ and B{X,T) similarly) in a sohition starting near 
the subspacc of non-modvilatcd solutions and converging to the periodic orbit. The 
top graph shows the time evolution of |^o| (dashed) and \Bq\ (solid). The other 
three graphs show \ Aj\ (dashed) and \Bj\ (solid) for j = 1, 2, 3 respectively.Note the 
initial decay when the travel time T is short, followed by growth of the perturbation 
as T increases. Parameter values are L = 30, /xi = 10~^, fi2 = 0.1, q = 10"^, and 
the coefficients are as in figure 1. 

lies inside the unit circle). Finally we consider only the minimum wavemimber 
imin = 2n/L since this is the mode with the highest growth rate for the points 
in the (//i, //2) plane that we consider. Hence the eigenvalue condition yields a 
relationship (implicitly) between T and q which can be compared with numer- 
ical simulations. This relationship simplifies in the asymptotic regime where 
q is small and T is large. 

Substituting for e± — exp[(ao ± co)T] in S and computing its eigenvalues we 
find 



m± 
where 



262e2coT + ^^2g4coT + ^2 ^ y 



2{hl + g2)e2coT 



(48) 
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It might appear that the dominant term for large T is the first term in the 
expression for y^, proportional to q^e^'^^, but this is incorrect. In fact, the 
dominant terms in the limit of small q and large T are those proportional to 
-2g4coT rpj^-g igg^f^g ^jjg asymptotic relationship 



as g — > 0. For two combinations of L and /i.2, figure 9 compares this asymptotic 
relationship with the exact implicit T-q relationship implied by (48) setting 
m+ = 1, and with the results of numerical integrations of the PDEs (6) - 
(7). In both cases there is clear agreement as long as q is sufficiently small. 
At small q the numerical simulations become more difficult, due to the inter- 
mittent nature of the dynamics. Given the several approximations involved in 
the analytic estimate of T(g) the results of figure 9 are enouraging. The ma- 
jor discrepancy in the analysis is the consideration of only one wavenumber 
£min = 27r/L in the estimate. Implicitly we assume that the 'most unstable' 
eigenfunction direction near P+ is aligned exactly with the single Fourier mode 



For larger q numerical simulations did not converge to a stable periodic orbit, 
but instead remained spatiotemporally disordered even after large integration 
times. Although solutions still often spend considerable amounts of time near 
the spatially-uniform state these events occur only intermittently. At much 
larger q, numerical integrations of the PDEs show that the cycle still plays a 
role in local organisation of the spatiotemporally complicated dynamics. The 
spatial structure of solutions often resembles a series of fronts between intervals 
of points which remain near either P+ or P_ for a time, and then switch rapidly 
to a neighbourhood of the other equilibrium. This is illustrated in figure 10, 
particularly around X — 90 and X — lAO (where the solution is close to P+), 
and at X = 120 where the solution is close to P_. For the parameter values of 
the figure, the pure mode equilibria have amplitude \B\ = 0.2. Nearby spatial 
locations are often still well-correlated in time, see figure 11. 



7 Stable travelling waves and complex spatiotemporal dynamics 

At large positive yUi and ^2 the TW states are restabilised after a Hopf bifurca- 
tion involving a branch of modulated waves generated by the global bifurcation 
in which the heteroclinic cycle loses stability. Surprisingly, for the spatially- 
extended system, these spatially-homogeneous TW states are also stable for 
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Fig. 9. The dependence on q of the approximate half-period T[q) of the periodic 
orbit near the heterochnic cycle, for two combinations of L and ^2- In both cases, 
the illustrative coefficient values were used, as for figure 1, and fii = q^- Squares 
(□) are the results of numerical simulations for L = 30 and /X2 = 0.1, and triangles 
(A) are numerical results for L = 25, ^2 = 0.05. Diamonds (O) give the exact T-q 
relationship implied by the linearised analysis near the cycle using I = imin = Stt/L, 
in both cases, and the asymptotic result of the linearised analysis is shown by the 
dashed and solid lines. 

large enough //i and /X2, as shown in figures 12 and 13. The stable TW coex- 
ist with stable complicated spatiotemporal dynamics, illustrated in figures 14 
and 15. The latter state might be expected to be a more generic solution in 
this region of the {^1,^2) plane, although no systematic study of the relative 
sizes of the basins of attraction has been performed. 
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Fig. 10. The moduli |^| (dashed hne) and \B\ (soHd line) for an instantaneous view 
of a typical solution of the PDEs (6) - (7). The horizontal dash-dotted line is at 
\B\ = \/\iJi2/o.2) = 0.2 corresponding to the amplitude of the P+ equilibrium. The 
parameter values are q = 0.2, = 0.04 and ^2 = 0.2; the illustrative coefficient 
values are used and the domain is of length L = 200 with periodic boundary condi- 
tions. 
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Fig. 11. Time evolution of at two spatial locations, X = 105 (solid line) and 
X = 110 (dashed line), for a typical solution of the PDEs (6) - (7). The parameter 
values, coefficients and domain size are as for figure 10. 
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Fig. 12. Temporal evolution of the spatial averages of the amplitude and phase 
variables at ni = 0.15, /X2 = 0.7, L = 30; all other coefficients are as in figure 1. The 
four panels show spatial averages of: (a) 1^41 (solid line) and \B\ (dashed line); (b) 
\Ax\ (solid) \Bx\ (dashed); (c) X (solid), 9 (dashed); (d) \Xx\ (solid), \9x\ (dashed). 
After an oscillatory transient the solution is attracted towards a spatially-periodic 
TW state where 9 (and cp, not shown) increase linearly with time while |^|, \B\ and 
X remain constant. 
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Fig. 13. The spatio-temporal evolution of the amplitudes |^| (upper panel) and \B 
(lower panel) for the parameter values of figure 12, showing convergence to TW. 
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Fig. 14. Temporal evolution of a persistent complex spatiotemporal state at 
Hi = 0.15, //2 = 0.7, L = 30; all other coefficients are as in figure 1. The four 

panels show spatial averages of: (a) \A\ (solid line) and \B\ (dashed line); (b) \Ax\ 
(sohd) \Bx\ (dashed); (c) X (solid), 9 (dashed); (d) \Xx\ (solid), \ex\ (dashed). This 
solution coexists stably with the TW state shown in figure 12. 
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time 

Fig. 15. Spatio-temporal evolution of the amplitudes \A\ (upper panel) and \B\ 
(lower panel) for the parameter values of figure 14, showing complicated spatiotem- 
poral dynamics. 
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8 Discussion and conclusions 



In this paper we have examined the robustness of the results of Part 1 on the 
dynamics near the 1 : 2 strong spatial resonance, to modulational instabilities 

which might be expected to play a role when the spatial resonance is close to, 
but not exactly, 1 : 2. We have derived a pair of coupled amplitude equations 
using symmetry arguments; these amplitude equations are PDEs rather than 
ODEs as is the case for exact resonance. 

We have examined the stability of the spatially-periodic solutions found in 
Part 1, looking for instabilities to finite- wavelength and long- wavelength per- 
turbations. In general, we find that the results of Part 1 still apply when ni 
and iJ,2 are large (equivalently when the mismatch q is small) , but that new in- 
stabilities are present near the codimension-two point i^i — H2 — 0. Since these 
lead, in general, to variations with a horizontal wavenumber 0{eq), solutions 
still resemble those of the non-modulated problem over short lengthscales. 

This paper is primarily concerned with identifying the several codimension- 
two points in the parameter space which link modulational instabilities and 
amplitude or phase instabilities which are present in the ODEs discussed in 
Part 1. We have presented reduced equations that describe the dynamics near 
these various points. The form of these reduced equations (which are still 
PDEs) may be deduced from symmetry arguments, and can also be derived 
through perturbation expansions. This analysis enables us to understand qual- 
itative features of the dynamics; although analytic solutions are still extremely 
difficult to derive for the reduced equations, this process of deriving 'normal 
forms' contributes greatly to an understanding of the interaction of these var- 
ious instabihtics. 

Our three main results are firstly, that the pure mode P and the mixed-mode 
M+ may undergo steady-state instabilities to finite-wavelength modes with a 
maximum critical wavenumber £c = q (for P) and £c = 2g (for M+). This 
confirms our intuition; it is natural to think that there are only two important 
lengthscales in the problem - the 0(1) lengthscale of the spatially-periodic 
patterns and long lengthscales > 0{l/eq) introduced by the deviation from 
exact resonance. M+ may also undergo an oscillatory bifurcation to a non-zero 
wavenumber mode, and numerical investigations suggest this wavenumber is 
also 0{q). Secondly, the reduction of the governing PDEs to a 'normal form' 
and subsequent analysis of these reduced equations for the codimension-two 
phase instability at the point X is novel and is presented in detail. Thirdly, we 
have considered the fate of the structurally stable heteroclinic cycle present in 
the ODE problem. Although the equilibria on the cycle, and hence the cycle 
itself, must be unstable to long-wavelength perturbations, the dynamics for 
small q result in a stable periodic orbit that hes close to the cycle and whose 
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period can be estimated well analytically. At larger q the heteroclinic cycle 
still organises the spatial dynamics over short distances, although solutions 
are spatiotemporally complex as might in general be expected. 

There is, of course, much more work that could be done on this problem. Two 
obvious directions of interest are the fate of the more complicated heteroclinic 
cycles examined by Porter & Knobloch [20] to modulational perturbations, 
and the behaviour for asymptotically large q (within this scaling ansatz); we 
have concentrated on behaviour for g < 1, where the results may be most 
easily interpreted with reference to previous work. 
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Appendix 



In this appendix we give the coefficients in the Ginzburg-Landau equation 
governing the instabihty of the pure mode solution P+ along the interior of 
the line OL, discussed in section 4. 



The P+ equilibrium is given by A = 0, B = y 1x2/0.2 = Bq. We start from the 
ansatz 



A{X, T) = e a{i, r)i?oe^^^ + a{i, r)2g(g + £)e 



-ilX 



7o(e, r) + 71 (e, r)e2^^^ + 72!^, r)e-^'« + 0{e'), (49) 



-2UX 



B{X, T)^Bo + [/?o(e, r) + M^, r)e'^" + T)e-''"J + 0{e'), (50) 

where ^ = eX is a new, longer, lengthscale and r = e^T is a new slower 
timescale. After a multiple-scales expansion, substituting this ansatz into the 
PDEs (6) - (7), the perturbation amplitudes Po, Pi, P2,lo,1i,l2 can be ex- 
pressed in terms of a. The usual solvability condition at O(e^) is satisfied 
identically, and at 0{e^) the solvability condition yields the Ginzburg-Landau 
amplitude equation for a: 



Ciar = Cia^^ — C2a\a\ , 
where the real coefficients Ci and C2 are 



62 = ai[B^ + 16q\q + if + IQBlq^q + if] 
+c^Bo[2WBl + Sq\q + ifh - 4g(g + £)] 
+Ah^Blq{q + £)(C4 + £5) - Blc^ - Aq\q + ifc^, 

C3 = -[B^oh + ^hq^q + if + Aq{q + £)]/(2a2Po), 

. ^ 2Blq{q + £)[a2g(g + ^) - 462c£^] - ^g(a2^g + 4c£^) 
""^^ 8c£2(a2Sg + 2c£2) 

. _ 02^0 - 2q{q + £)[g(g + f)(a2gg + cf) + 462Pgc£^] 
Sd^{a2Bl + 2cP) 

Hence C2 > when ai is sufficiently positive. In the special case 03 = 64 = 05 = 
0, the expression for 62 agrees with that derived by CouUet et al. [6]. In fact, the 
O(e^) amplitudes 70,71,72 in (49) do not contribute to the complexity of the 
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expression for C2; the complexity is derived from the third-order contributions 
from the B equation. 
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